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I. INTRODUCTION 



The concept of symmetry has played a central role in modern theoretical physics. A 
discovery of a new symmetry enlarges our ability and possibility to describe new phenom- 
ena both in the physical nature and in mathematical models. Conversely, it may be almost 
certain that there is an underlying symmetry if a system under consideration exhibits a 
significant property that is not shared in the generic cases. Actually, A/'-fold supersymme- 
try was discovered from the observation of the disappearance of the leading divergence of 
the perturbation series for the specific energy levels of a quantum mechanical model at the 
particular values of a parameter involved in the model This symmetry is ageneraliza- 
tion of the ordinary supersymmetry in one- dimensional quantum mechanics [2, 01 and is 
characterized by A/'th order derivative supercharges. Similar gen eralizations were found and 
investigated in various different contexts 0, II S Q i, i, HO, 111 H Q Q El [B, H Q ■ 
A generalization to A/'th order derivative supercharges for all M is, however, difficult and 
most of the investigations by the other authors were limited to the cases of the second-order 
generalization and to the cases constructed by the factorization method. Recently, we found 
another A/'-fold supersymmetric model ^3 ^'^^ further found general forms of an AA-fold 
supersymmetric family which we have called type A j23]. On the other hand, several quasi- 
solvable models j2l|, 22| were found to have A/'-fold supersymmetry in Refs. jH E3. 125I . l26l | 
in addition to those reported in Refs. 0,0]- Then, it was proved in generic wav |27l| that 
A/'-fold supersymmetry is essentially equivalent to quasi-solvability. Furthermore, it was also 
proved j28| that the type A A/'-fold supersymmetric models are essentially equivalent to the 
quasi-solvable models constructed by st(2) generators ji^. Dynamical properties of the A/'- 
fold supersymmetric models were discussed in Ref. ^3 ^^^1 investigated for a couple of mod- 
els in Ref. |31| . The nonperturbative analyses carried out in Ref. ^] together with those in 
Ref. 0] revealed several significant properties that the type A A/'-fold supersymmetric mod- 
els share. Especially, we clarified the important role of the normalizability of the solvable 
sector, which is crucial for the dynamical A/'-fold supersymmetry breaking but has rarely dis- 
cussed by the other authors. Furthermore, it was argued that A/'-fold supersymmetry is not 
only sufficient but may also be necessary for the existence of convergent perturbation series. 
Up to now, most of the quasi-solvable one- dimensional quantum mechanical systems belong 
to the sl(2) quasi-solvable models. Several new findings have been reported in connection 
with them. One of the examples is the new orthogonal polynomials firstly found by Bender 
and Dunne [s^- The idea was soon generalized to all the sl(2) quasi-exactly solvable mod- 
els (33I Is^ l . Furthermore, several realistic physical systems have been found, which can be 
reduced to one-dimensional quasi-solvable models 0, IH IH 113, HI H il 41, il il 4^ . 

In this article, we will report the general aspects of type A A/'-fold supersymmetry, most 
of which remain unsolved yet. In Ref. it was discussed that, in view of the sl(2) con- 
struction of type A A/'-fold supersymmetry, the condition for type A A/'-fold supersymmetry 
derived in Ref. |0| provides only a sufficient but not necessary conditions. On the other 
hand, only one of the pair of the type A Hamiltonians was investigated in Ref. ji^ and 
relations between the pair in view of the 51(2) structure have not been known. In other 
words, the equivalence between the analytic construction in Ref. [23] and the algebraic one 
in Ref. [2^ has not been established. Another problem to be solved is the classification 
problem. In Ref. [s^l, we attempted to classify the type A A/'-fold supersymmetric models 
by considering the invariance under the linear transformations but found to be incomplete. 
Later we found that the complete classification of the s[(2) quasi-solvable models was already 
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achieved in Ref. |45|,|46| by the consideration of more extensive GL{2, R) transformations. 
Therefore, we have recognized that the complete classification of the type A models should 
be done by investigating the invariance of both the Hamiltonians and A/'-fold supercharge 
under the GL(2,M) transformations. Lastly, the structure of the anti-commutator of the 
type A A/'-fold supercharges, which we have called mother Hamiltonian, has not been inves- 
tigated yet at all and thus we have not established the complete type A A/"-fold superalgebra. 
We will give the complete answers to the above problems in this article. 

The article is organized as follows. In the next section, we give a brief review on the 
definition and the general properties of A/'-fold supersymmetry needed in this article. In 
Section IIIIl we define type A A/'-fold supersymmetry and give two different approaches, 
namely, analytic and algebraic approaches to construct the type A A/'-fold supersymmetric 
models. The classification problem is discussed in Section HVl By considering the invariance 
property of both the Hamiltonians and A/'-fold supercharge, we obtain the complete answer 
to the problem. In Section |3 we investigate 2- fold supersymmetry with an emphasis on the 
uniqueness of supercharges and the polynomiality of the 2-fold superalgebras. A novel fea- 
ture of weak quasi- solvability defined in Section |TT1 is discussed. In Section IVTl we first show 
the polynomiality of the type A A/'-fold superalgebras. Then, by a suitable generalization 
of the original idea of the Bender-Dunne polynomials to all the type A A/'-fold supersym- 
metric models, we show that the anti-commutators of the type A A/'-fold supercharges are 
expressed as the critical polynomials of the generalized Bender-Dunne polynomials. We 
claim that both the normalizability of the solvable sector and the weak orthogonality of the 
polynomials do not play essential roles on the generalization. A novel interpretation of the 
critical polynomials is given in view of polynomial invariants. Interesting future problems 
are discussed in the final section. In Appendix we summarize the results of the invariant 
theory of polynomial systems needed in this article. 

II. A/'-FOLD SUPERSYMMETRY 
A. Definition 

We start with defining A/'-fold supersymmetry in one dimensional quantum mechanics. 
Let us introduce a bosonic coordinate q and fermionic coordinates ip and ip"^ satisfying. 



(2.1) 



Hamiltonian H^v" is given by. 



(2.2) 



where are ordinary scalar Hamiltonians: 



(2.3) 



with p = —id/dq. A/'-fold supercharges Qj^ and Qj^- are introduced by. 



(2.4) 
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where P^/" is given by a polynomial of TVth degree in p: 

Pn = P^ + WM-i{q)p^~^ + ■■■ + wi{q)p + Wo{q), (2.5) 

that is, Pj\f is an A/'th-order linear differential operator. Then, the system (j2.2j) is defined 
to be M-fold supersymmetric if the following algebra holds: 

{Q^,QAr} = {gir>QV} = 0, (2.6) 
[Q^,H^] = [gV,H^^] = 0. (2.7) 

The former relation is trivial due to Eq. ()2.H) while the latter is equivalent to the following 
intertwining relations: 

PmH]^ - H+^Pm = 0, Pj^H^ - H^Pl^ = 0. (2.8) 

Therefore, the relation (j2.8p gives the condition for the system to be A^-fold supersym- 
metric. From the definition, it is evident that A/'-fold supersymmetry reduces to the ordinary 
supersymmetry HQ when A/" = 1. 



B. Quasi-solvability 

The A/'-fold supersymmetric models defined above have several significant properties sim- 
ilar to those of the ordinary supersymmetric models [l^]. One of the most notable ones is 
quasi-solvability j^, i^l- A linear differential operator H of a single variable q is said to be 



quasi-solvable if it preserves a finite dimensional functional space VV whose basis admits an 
explicit analytic form: 

HVm C Vat, dim Vat = n{M) < oo, = span {0i(g), . . . , 0„(Ar)(g)} • (2.9) 

An immediate consequence of the above definition of quasi-solvability is that, since we can 
calculate finite dimensional matrix elements 8^,/ defined by, 

n{Af) 

= 5^ Sfc,,0; {k = l,...,n{X)), (2.10) 
(=1 

we can diagonalize the operator H and obtain the spectra of it in the space VV, at least, 
algebraically. Furthermore, if the space VV is a subspace of a Hilbert space L^(S') {S C M) 
on which the operator H is naturally defined, the solvable spectra and the corresponding 
vectors of VV give the exact eigenvalues and eigenfunctions of H, respectively. In this case, 
the operator H is said to be quasi- exactly solvable. The role of the normalizability of the 
solvable sector is investigated in view of dynamical properties in Ref . 0, |^ . To construct a 
quasi-solvable model, it is convenient to introduce an A/'th-order linear differential operator 
P and define the vector space as, 

V^ = kerP. (2.11) 

Now, it is easy to see that an operator H is quasi-solvable with the solvable sector ()2.1H) if 
the following quasi-solvability condition holds 2ol |: 



PHVm = 0. (2.12) 
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This formulation enables us to clarify the relation between quasi- solvability and A^-fold 
super symmetry. Indeed, it can be easily shown that all the quasi-solvable models satisfying 
Eq. ()2.12|) are A/"- fold supersymmetric if we set = P, Hjf = H and Hj^ = H + iw'j^_^, 
where wj^^i is defined in Eq. (12. 5p . The converse is also true. From the intertwining 
relations (j2.8p . we find that all the A/'-fold supersymmetric systems are quasi-solvable: the 
quasi-solvability condition ()2.12|) holds for H = H^{H^) and P = Pj\f{Plf), respectively. 
Therefore, if we define, 

Vjj- = her Pj^ = span : n = 1, . . . ,M}, (2.13a) 
= kerP| = span{0+ : n = 1, . . . ,Af}, (2.13b) 

we have H^V^ C and the following Schrodinger equations on the subspaces V^: 

M 

hM = E ^ti^^ = ^t4>t (A; = 1, . . . , AT). (2.14) 

1=1 

The spectra are determined from the characteristic equations for the matrices S^: 

detM^(P±) = 0, M^(A) = 2(AI- S^). (2.15) 

We should note that, for a given operator P, we cannot always obtain analytic solutions 
of Eq. ()2.11|) . Therefore, quasi-solvability formulated from Eqs. ()2.11|) and ()2.12j) is less 
restrictive than the one defined by Eq. (j2.9|) . In the situation where this difference is crucial, 
we may be better to call the less restrictive case weakly quasi-solvable. With this terminology, 
we say more correctly that A/'-fold supersymmetry is equivalent to weak quasi-solvability. In 
Section El we will discuss weak quasi-solvability again. 



C. Mother Hamiltonian 

In the ordinary supersymmetry, the anti-commutator of the supercharges corresponds to 
the Hamiltonian. However, it is not the case in A/'-fold supersymmetry. This is because 
{Qjy, Qa^} is now a 2A/'th-order differential operator. The half of the anti-commutator is 
called mother Hamiltonian and is denoted by 'H^f. 

'Hn = \{Q^m^Qn}- (2.16) 

An immediate consequence of the above definition is that the mother Hamiltonian always 
commutes with the A/'-fold supercharges, that is, it is A/'-fold supersymmetric: 

[QAr,^Ar] = [gJv-,^Ar] = 0. (2.17) 

Furthermore, if the original Hamiltonian Hyv" is A/'-fold supersymmetric, the mother Hamil- 
tonian also commutes with H_/v^ due to the relation ()2.7|1 : 

[HAr,HA^] = 0. (2.18) 

From the above relations, it is expected that the mother Hamiltonian Tij^ has an intimate 
relation with the original Hamiltonian H_^/-. Indeed, it was shown j2^ that, if the A/'-fold 
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supercharges Qj^ and Qjy- are uniquely determined and Hj^ ^ Hj^, H^f is expressed as the 
characteristic polynomial of A/'-th degree for appeared in Eq. ()2.15|) with the argument 
replaced by H_^: 

nj^ = l det Mj^iU^f) = ^ det M^(H^). (2.19) 

III. TYPE A TV-FOLD SUPERSYMMETRIC MODELS 

In contrast to the ordinary supersymmetric quantum mechanics, the construction of an 
A/'-fold supersymmetric model is a non-trivial problem. In the case of ordinary supersym- 
metry, the mother Hamiltonian ()2.16|) defined by the anti-commutator of the supercharges 
is a desirable Schrodinger operator of the form (j2.2p and thus can be identified as a super- 
symmetric Hamiltonian due to the relation (j2.17j) . In the case of 2- fold supersymmetry, the 
condition (|2.8|) can be completely solved and thus the most general form of the 2-fold su- 
persymmetric Hamiltonian and the 2- fold supercharge have been known 0,123]. However, 
we can hardly solve the condition ()2.8|) for A/" > 3 in general. Nevertheless, we have found 
a special case where the condition ()2.8p can be solved for arbitrary A/". We have called this 
case type A [20] . The type A A/'-fold supercharge is defined as the following special form of 
A/'th-order linear differential operator: 

Pm= n {p-iW{q)+ikE{q)) 

k=-{M-l)/2 

p - iW{q) + i '^ ^ ^ E{q)^ (^p - iW{q) + i '^ ^ ^ E{q)^ x ■ ■ • 
■X (p- iW{q) - (p - iW{q) - z^^^^E(g)^ . (3.1) 



In the above definition, we note that E(q) and W{q) correspond to E{q) and W{q), re- 
spectively, in the previous articles ji^ 113, lisL l30| . We will use tildes for another particular 
purpose (see Eqs. ()3.3|) and below) and thus we do not follow the old notation anymore. A 
system fl2.2|) is said to be type A M-fold supersymmetric if the condition ()2.8|) is fulfilled with 
this type A A/'-fold supercharge. There are two ways to construct a type A model, namely, 
analytic and algebraic constructions. The former is to solve the intertwining relation (j2.8p 
directly while the latter is to solve the quasi-solvability condition (j2.12j) . In the following 
sections, we will first show the analytic construction and next the algebraic one. 



A. Analytic Construction 



The condition for type A A/'-fold supersymmetry was firstly investigated with the aid of 
induction in Ref . Il^ . Later, it was reexamined in the context of the algebraic construction 

gives a 



isl ). Then, it has turned out that the set of the conditions derived in Ref. 2^1 only 



sufficient one. The origin of the defect was explained in Ref. 2^]. In the following, we will 
give an improved direct proof of the necessary and sufficient condition for type A A/'-fold 
supersymmetry. 
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The necessary and sufficient condition for the system ()2.2j) to be A^-fold supersymmetry 
with respect to the type A A/'-fold supercharge ()3.1|) is the following: 

V^{q) = \w{qf + \v^{q\ v^{q) = {E{qf - 2E'{q)) ± UW\q) - 2R, (3.2a) 

- ^(^)) ^ + ^^^)) ^(^) = ^ 2, (3.2b) 

(A-2i.(,))(^-i.(,))^(^ + E(,))E(,)^0 for Ar>3, (3.2c) 
where R is an arbitrary constant. 

Proof. The proposition will be proved by induction. At ffist, we make gauge transformations 
on Pa/ and to facilitate the calculations, as follows: 

iV = t^{GxU)PM{G^U)-^ = n " ^^)' (3.3a) 

fc=0 

= {GmU)H^{GmU)-\ (3.3b) 

where Gj^ and f/ are defined by, 

Gm = exp (^^^^ JdqE{q)^ , f/ = exp (^jdqW{q)^ . (3.4) 

In the above and hereafter, we attach tildes to operators, vectors and vector spaces to 
indicate that they are quantities gauge-transformed with the gauge factor Gj^U . We set, 

I^ = 2{P^H^-H+P^). (3.5) 

We assume that = for a natural number Af if the set of the conditions (j3.2p holds for 
this J\f. Then, we will prove that /at+i = if and only if the set of the conditions (j3.2|) holds 
for A/" replaced by A/" + 1. From the following relation: 

2{Gm+iU)H^{Gm+iU)-' = 2H^ + Ed+ ^E' - ^:^^E^ - EW, (3.6) 



we obtain, 



2H^^, = {Gu+iU){2H^ -v^ + v^^,){Gu+iU)-^ 

1 ^, 2U-1 
2 

2H^ + Ed + u%^,-u%, (3.7) 



2H^ + Ed + -E' -—E'-EW-v^ + v^^, 



where are defined by. 



u 



, .AT-l , (AT-l) 



2 



,^ = v^ + W' + -^^E' - ^ -{AT- l)EW. (3.8) 
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In the above, we note that are given by Eq. ()3.2a|) from the inductive assumption while 
v^^i are unknown functions to be determined. From Eq. ()3.7p . we obtain, 

Ij^+i = Pn+1 (2^^ + Ed + uj^^, - M^) - (2H+ +Ed + ^1 - iV+i 

+ Pu+iEd - EdpM+i + [iV+i, - t^jv-]- (3-9) 
The last four terms in the r.h.s. of Eq. ()3.9p are calculated as follows: 

2[H^, d - ME] = ((AT + l)E' - 2W') (d - UE) 

- {uXf + 2UWE -ME"^ -NE')' , (3.10a) 

A/'+l N 

Pu+iEd = EP^+2 + kE)E^.^)P^, 

n=l k=n 

MM 

= ePm+2 + (AT + + XI ^ n kE)E^o,-i)Pn, 

n=l k=n+l 

(3.10b) 

EdP^^i = EP^+2 + i^f+ 1)E^Pm+i, (3.10c) 
M M 

[Pm+i,uJ^+i -ujf\ = ^ W {d- kE){u^^-^ - M^)(o)Pn, (3.10d) 

n=0 k=n+l 

In the above and hereafter, we employ the following abbreviations: 

/(fc) = id-kE)f, /(fc,...) = (9- A;E)/(..,). (3.11) 

The following formula is useful for the calculations: 

M M-l M~l M n'-l 

l[{d-kE)f = fl[{d-kE) + Y, n id-kE)f^i^l[{d-k'E). (3.12) 

k=n+l k=n n'=nk=n'+l+l k'=n 

Substituting Eqs. ()3.10p for Eq. (j3.9p . we have, 

J^+i = {2W' + M^+i -ulf- u+^^^ + u^) Pm+1 + {u^ + 2AriyE - ME'' - ME') ' P^ 

MM 

+ E n - [("A^+l - V)(0) + ^^(0,-1)] Pn- (3.13) 

n=0 k=n+l 

From Eq. ()3.13p . we see that /v+i is at most (7V+l)th-order differential operator. Therefore, 
Jv+i = if and only if all the coefficients oi (fc = 0, 1, . . . , A/" + 1) vanish. Only the first 
term in the r.h.s. of Eq. ()3.13|1 contains the d^^^ term. Therefore, one of the conditions for 

J^+i = reads, 

2W' + M^+i - - M^+i + = 0. (3. 14) 

Combining Eq. ()3.14|) with Eqs. ()3.2a|l and ()3.8|1 . we have, 

v^^,-v-^^,=2{M+l)W'. (3.15) 
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Applying the formula ()3.12p and the condition ()3.14p to Eq. ()3.13p . we obtain, 
'AfiAT+l) 



-E(o,-i) + (AT + l){uj^^, - Mj^)(o) + {uj^ + 2MWE - ME^ - ME') 



P. 



M 



n (-9- fc^) [(^^^+1 -w^) (1,0) +^^(1,0,-1)] ^n'. 

n=0 n'=n fc=n'+2 



(3.16) 



Only the first term in the r.h.s. of Eq. (j3.16p contains the term. Therefore, we obtain 
another condition for /at+i = 0: 



E(^o,^i) + (AT + - V)(o) + {um + 2Ariy^ - ^fE' - ME')' = 0. (3.17) 



Combining the assumption ()3.2a|) with Eqs. ()3.15|) and ()3.17|) . we finally have, 
± (Ar+ 1)2-1 



12 



{E^ - 2E') ± (TV + 1)W' + const. 



(3.18) 



The resulting f^+i is nothing but the assumed form ()3.2aj] with A/" replaced hy J\f + 1. 
Therefore, the first condition ()3.2a|) has been proved inductively. Under the above conditions 
satisfied, we have from Eqs. ()3.8p and ()3.18|) . 
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(3.19) 



Substituting the above for Eq. ()3.16|) . we obtain. 



w=-EE n (^-^^) 

n=0 n'=n k=n'+2 



^(1,0,-1) - 



3n - AA + 1 



-^(1,0,-1) 



Pn'- 



(3.20) 



The above expression for is further arranged with the aid of the formula ()3.12|1 as 

follows: 



AT M-1 



n— n n' — n ^ 



3n - AA + 1 



n=0 n'=n 



n=0 n'=n n"=n' fc=?i"+3 

N{N + 1) 



-^'(i.o -1) J -Pa/^-i 

-EEE n (a-*^)(wWo,-«-^^^^fii 

W^(i,o -i)-fV-i 



-^(2,1,0,-1) -P» 



2 

^ A/'-2 M , 

- 9 + + 2) n - (^(2,1 A-i) - 



n - A/'+ 1 



-^(2,1,0,-1) -Pn- 



n=0 



fc=n+3 



(3.21) 
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When AT = 1, J2 = — Vr(i,o,-i) and thus we obtain, 

/. = « w,,,..,^{^^-E)±(^±^Eyv = 0, (3.22) 

i.e., the condition ()3.2b|l in addition to Eq. ()3.2a|) . When Af > 2, the condition ()3.2b|) has 
been aheady assumed and thus. 



^ A/-2 A/ 

W = g$^(^ + l)(n + 2) H {d - kE)E^2,i,o,-i)Pn. (3.23) 

n=0 fe=n+3 

When Af = 2, I3 = — -£"(2, 1,0,-1) /3 and thus we obtain, 

i.e., the condition ()3.2cj) in addition to Eqs. ()3.2a|l and ()3.2b|l . When J\f > 3, the condition 
fl3.2cj) has been aheady assumed too and thus, 

Im+i = 0. (3.25) 
Therefore, no additional condition is required any more. □ 

From the set of the conditions ()3.2|) . we find a procedure to construct a type A A/"-fold 
supersymmetric model as follows: 

1. Find out a particular solution E{q) of the nonlinear differential equation ()3.2c|) . 

2. Substitute the above E{q) for the general solution of the linear differential equation 
()3.2b|) given by, 

iy(g) = Cie-/'^''^('?) Jdq (^e/'^''^^'?) Jdqel"'^''^'^^^ 

+ C2 e- ^'^"^('^^ J dq e^'^^^'^^ + C3 e' ^'^^^^'^^ {d = arbitrary constants) , 

(3.26) 

to obtain W{q). 

3. Substitute the above E{q) and W{q) for Eq. ()3.2a|) to obtain the pair of potentials 



B. s[(2) Construction 

As previously noted, AT-fold supersymmetry is essentially equivalent to quasi-solvability. 
Recently, some special quasi-solvable models which can be constructed from the sl(2) gen- 
erators jlil were found to be type A A/"- fold supersymmetric [l[ 113, [i^ [i^, 0, . 
Then, it was shown in Ref. [2^ that the type A A/'-fold supersymmetric model is essentially 
equivalent to the s[(2) quasi-solvable model. In this section, we will review the equivalence 
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and then clarify the relation between the pair of the Hamiltonian in the framework of 
the sl(2) quasi-solvable models, which has not been discussed yet in the previous articles. 
Let us first construct Hjj- so that it is quasi-solvable with respect to the type A operator 
(|3.1|) . The quasi-solvability condition (|2.12p in this case is, 

PmH]^V]^ = 0, (3.27) 

where Vj^ is given by Eq. (|2.13a|) . On the gauge-transformed space, the condition (|3.27|) is 
equivalent to, 

PmH^Vj^ = 0, (3.28) 
where Pj^ and Hj^ are given by Eqs. (j3.3|) and is defined by, 

Vj^ = ker Pj^ = span {Gj^U(f) : G ker Pj^} . (3.29) 
Introducing a function h{q) defined as a solution of the following differential equation: 

h"{q)-E{q)h'{q)=0, (3.30) 
we find that Pj\f is expressed in terms of h as, 

/'^-Or^. (3.31) 

Thus, we easily have, 

= span {1, h{q), . . . , /i(g)^-^} . (3.32) 
From Eqs. ()3.28|) and ()3.32|) . the quasi-solvability condition for Hj^ reads. 



i7-/ife-i = for V/c = 1,...,A/'. (3.33) 



Any constant is a trivial solution of (|3.33p . For TV = 1, any first-order differential operator 
of the following form, 

fm^, (3.34) 

is a solution, while for A/" > 2 there are three independent first-order differential operators 
as solutions of (j3.33j) : 

The J+'°'^ defined above satisfy the s[(2) algebra: 

[J+, J-] = -2J°, [J±, J°] = tJ^. (3.36) 
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In the same way, we find that for = 1,2, any second-order differential operator of the 
following form. 



f2{h) 



(3.37) 



is a solution, while for A/" > 3 there are five independent second-order differential operators 
as solutions of ()3.33p : 



h 



. d^ 



I T-\2 , TO T- A/" - 1 _ 

(^) > ^^ = ^°^ 

, (AA-l)(AA-3) 



(J°)' + (Ar-2)J° + 



/.sil _ (AT - l)(Ar - 2)/i = J+ + J+, 



dh^ 



1{M -2)h^^^{M -\){M -2)h^ = (J+y 
dh ' 



Therefore, the general solution of ()3.33p for A/" > 3 can be expressed as. 



H 



(3.38) 



where a-^ , b\ and •* are constants. This gauged Hamiltonian is nothing but the sl(2) 

quasi-solvable model [i^. We can set = without any loss of generality due to the 
relation: 



1 



(J+ J- + J- J+) - (J' 



0\2 



1 



2' ' ^ ' ^ 

For convenience, we introduce new parameters as follows: 



r{-) 



6i = -6^-\ 60 



,(-) 



Then, the gauged Hamiltonian ()3.38|1 is expressed in terms of h as, 

d^ 



(3.39) 

(3.40a) 
(3.40b) 



R 



N 



dK^ 



Af -2 
Qih) - -^^P'ih) 



d_ 

dh 
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where P{h) and Q{h) are given by. 



Q(/i) = b2h^ + bih + bo, 



(3.41) 



(3.42) 
(3.43) 
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while R = C*^~^ + (A/'^ — 1)02/12 is a constant. If the Hamiltonian (j3.4H) is gauge-transformed 
back to the original Hamiltonian with the gauge factor G^U, it is in general not of the 
canonical form of the Schrodinger operator like Eq. ()2.3j) . We can find that the original 
Hamiltonian Hj^ becomes of the canonical form if and only if the following conditions hold: 

P{h) = l{h\q))\ (3.44a) 
Q{h) = -W{q)h'{q). (3.44b) 
Under the above conditions satisfied, we have, 

= {G^U)-^H^{GmU) = -^-^ + Vj^iq), (3.45) 



where the potential Vj^ is given by, 
1 



UP 



{N^ - 1) (^PP" - \{P'?^ + 3J\f{P'Q - 2PQ') - 3Q2 



- R. (3.46) 



Substituting Eqs. ()3.3Up and ()3.44|) for Eq. ()3.46|) . we can check that the above expression 
fl3.4fi|l for Vj^ is identical with the one in Eq. ()3.2aj) . 

When M = 2, any second-order operator ()3.37j) can be added to the gauged Hamiltonian 
(I3.38|l and thus all the Uij in Eq. ()3.38j) can be arbitrary functions of h. As a consequence, 
P{h) in Eq. 1)3.411) can be an arbitrary function of h. 

When jV = 1, any first-order operator ()3.34p in addition to any second-order one ()3.37|) 
can be added to the gauged Hamiltonian ()3.38p and thus all the bi in addition to all the a^j 
in Eq. ()3.38p can be arbitrary functions of h. As a consequence, both P{h) and Q{h) in 
Eq. ()3.41|) can be arbitrary functions of h. 

Combining the above considerations for TV = 1,2 with the results ()3.42j) and ()3.43j) for 
A/" > 3, we obtain the following conditions for P{h) and Q{h): 

-^Q{h)=0 for Af>2, (3.47) 

—rPih) = for ^^ > 3. (3.48) 
dh° 

If we rewrite the above conditions in terms of q with the aid of Eqs. ()3.30)1 and ()3.44j) . 
we find that the condition (|3.47p is equivalent to Eq. ()3.2b|) while the condition ()3.48j) is 
equivalent to Eq. (|3.2cp . Therefore, the complete correspondence between type A A/"- fold 
supersymmetric models and the sl{2) quasi-solvable models has been established except for 
the other Hamiltonian Hj^, which will be discussed in the next. 

The partner Hamiltonian Hj^ should be constructed such that it is quasi-solvable with 

jy. 11 IS cixay lu sec thai iiic upciaiui l 



respect to Ptr. It is easy to see that the operator Pjr can be converted into the form identical 



-1. 



with Pj^ by a gauge transformation with another gauge factor Gj^U 

Plf = {GmU-') P| {G^U-r' = {h'f (3.49) 

where h is the same function of q as defined previously by Eq. ()3.30|) . In the above and 
hereafter, we attach bars to operators, vectors and vector spaces to indicate that they are 



13 



quantities gauge-transformed with the gauge factor Gj^U~^. Therefore, the quasi-solvabihty 
condition for the gauged Hamiltonian is completely the same as that for Hjj-, Eq. ()3.33p . 
This means that is another quasi-solvable model constructed from the s[(2) generators 
()3.35p and thus has the same form as Eqs. (j3.38|) and (j3.41|) : 



J2 a^^\j\P+ J2 t^J'-C^^^ 



(3.50a) 



-P^h) 



i=+,0,- 



Af 



P+\h) 



d_ 

dh 



(3.50b) 



where a'^\ b^^\ C(+) and R+ are constants, and P^{h) and Q~^{h) is a polynomial of fourth- 



ly ' I ' 



and second- degree, respectively. Up to now, there is no relation between and Hjf. If we 
demand that be gauge-transformed, with the gauge factor Gj\/U~^, back to an operator 
of the Schrodinger form, we have. 



P+{h) = -{h\q)f, Q+ih)=Wiq)h'iq). 



Combining the above with Eq. ()3.44|1 . we yield the following relations: 

P+{h)=P{h), Q+{h) = -Q{h). 
Then, the partner Hamiltonian Hj^ becomes. 



2dq^ 



+ V^{q), 



with, 



UP 



(M^ - 1) ( PP" - ^{P 



t\2 



3M{P'Q - 2PQ') - 3g2 



(3.51) 



(3.52) 



(3.53) 



(3.54) 



Af 



To establish the relation between R and R^, we should recall the formula j2j 
iw'j^_^, where wj^^i is defined in Eq. (|2.5|) . For the type A A/"- fold supercharge (j3.1|) . we 
have wj^^i = -iJ\fW. From Eqs. (I3.44bjl . and ^^)^ we finally obtain. 



R+ = R. 



(3.55) 



We can easily see that the potential ()3.54|) with Eq. ()3.55|) is identical with in Eq. ()3.2ap . 
Thus, the models constructed in this section are completely identical with the ones by the 
analytic construction in the preceding section. 

Before closing the section, we will refer to an interesting relation between the algebraic 
Hamiltonians Hjj- and Hj^. From Eqs. ()3.52|) and ()3.55p . the relation between the parameters 



in Hj^ and the ones in reads. 



+ bi^ + a 



(-) _ 
+0 - 



,(-) 



(Vz,j = +,0,-), C(+) = C(-), 



0, 



^0- 



0. 



(3.56a) 
(3.56b) 
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Substituting the above for Eq. ()3.50a|) . we can rewrite in terms of a-^ and h\ . 

H^ = - ^^J'J'- E t^J'-C^'^- (3.57) 

jj=+,0- i=+,0,- 

Comparing this expression with Eq. ()3.38|) . we see that H'^ is obtained from Hjj- by inter- 
changing the order of the quadratic terms and by interchanging the sign of the coefficients 
of the hnear terms. 



IV. CLASSIFICATION OF TYPE A MODELS 



It was shown that the s[{2) quasi-solvable models can be classified using the shape in- 
variance of the Hamiltonian under the action of GL{2, K) {K = M or C) of linear fractional 
transformations 4^ 4^. The equivalence established in the previous section ensures that the 



type A AT-fold supersymmetric Hamiltonians can be fit into the same classification scheme. 
Then, a natural question arises whether the type A A/"-fold supercharges can be also classi- 
fied in the same scheme or not. This question was left as an open problem in the previous 
paper [s^ and will be completely answered in this section. 



A. GL{2, K) Shape Invariance 

The linear fractional transformation of h is introduced by, 

h^h = ?^h±^ /?, 7, 5 G iT; A = a(5 - /37 ^ 0). (4.1) 

Then, it turns out that the gauged Hamiltonian ()3.4H1 is shape invariant under the following 
transformation induced by Eq. ()4.1|) : 

H-^{h) ^ Hj^ih) = i^h + 5f-'H^ih)ilh + 5)-^^-'\ (4.2) 
where P{h) and Q{h) in Eq. ()3.41|) are transformed according to, 

P{h) ^ P{h) = A-\-fh + 6)^P{h), (4.3a) 
Q{h)^Q{h) = A-\^h + 6fQCh). (4.3b) 

B. Invariance of the Hamiltonian 

In the next, we will examine the transformation of the original Hamiltonian H^. Since 
the function h{q) is determined by Eq. (j3.44aj) . we have, 

, , [ dh , , f dh /"„ A 1 

\q\ = / — 1-^ |g| = / — = dh 



dh 



2P{h) 



(4.4) 
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Therefore, h{q) before the transformation is identical with h{q) as a function of q. In other 
words, the relation between h before the transformation (denoted by hold) and h after the 
transformation (denoted by knew) is, as a function of q, consistent with Eq. ()4.1|1 : 

lKew[q) + 

On the other hand, the potentials can be rewritten as j4 

= [2{U' - l)H[P] T Qf^ - QQ'] - R, (4.6) 

where H[P] is the Hessian of P and {P,Q)^^^ are the first transvectant of P and Q, given 
by (see also Appendix . 

H{P) = PP" - ^ {Pf , (P, Q)(^) = 2P'Q - 4PQ'. (4.7) 

All the objects 0{h) in both the numerator and denominator of Eq. ()4.f)|l . namely, H[P], 
{P,Q)^^\ and P transform according to, 

OiKid) ^ d{K,^) = A-'i^K,^ + 6yoCh), (4.8) 

that is, they belong to the multiplier representation p4, -2 of GL{2,K) defined in Eq. (jAlj) . 
As a consequence, the functional form of are preserved under the transformation: 

Finally, since the functional form of hoid{q) and h{q) is identical with each other, Eq. ()4.5p . 



the potential V)^ is invariant as functions of q 



V^ihoiM) ^ V^ihneM) = V^{h{q)) = V^iKidiq)). (4.10) 

Therefore, the type A A/'-fold supersymmetric Hamiltonians are invariant under the 
GL{2, K) transformation. 



C. Invariance of the Supercharge 



In the next, we will examine the transformation of the type A A/'-fold supercharge Pj^. 
The gauged A^-fold supercharge (|3.3ip is transformed according to. 



{h 



On the other hand, the gauge factors Gj\fU'^^ are transformed as. 



(4.11) 



GxU^^{hoid) = exp 



dq 



M-l 



E±W 



exp 



dh 



old' 



{M-l)P'{Kid)T'2Q{k 



old 



GmU^\K 



exp 



dhr, 



^P{hoid 

{M -l)P'{hne^)T'2Q{hnev^) 



(4.12) 
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From Eqs. (P?T|) and (jOj), the r.h.s. of Eq. (jil^ is calculated as 

Gj^U^^{hnew) = exp 



AP{h) J iKew + 

= {iKe^ + 5f-'GMU^\h). (4.13) 
The A/'-fold supercharge Pa/ is expressed as, 

Pn = {-i)^{GMU{Kid)Y^PN{GMU{hM)) 

= [-^h',^,f{GMU{hM)) -^{GMU{Ku)). (4.14) 

Combining Eqs. ()4.11|) - ()4.13|) . we obtain, 

Pm^Pm= Hf{G^{hne^)y'PM{G^{hne^)) 

= {-^h'f{G^fU{h))-'^{GuU{h)). (4.15) 

Comparing Eq. ()4.14|) with ()4.15|) . we see that Pj^ is obtained from Pj^ with hoid{q) replaced 
by h{q). From the fact that h{q) is identical with hoid{q) as a function of g, Eq. ()4.5|1 . we 
finally conclude, 

Pn = Pm, (4.16) 

that is, the A^-fold supercharge is also invariant under the GL{2, K) transformation. We 
note that the manifest invariance of the A/'-fold supercharge is lost if we express it in terms 
of E{q) and W{q) as Eq. ()3.H) . This is because E{q) is not an invariant function under the 
GL{2,K) transformation, as we will see below. From Eq. ()4.5|) . we have, 

Kidiq) ^ ^ h';^eM _ '^iKeM ^ 17) 

The function h{q) is defined so that Eq. ()3.30p is fulfilled and thus the relation between E{q) 
before and after the transformation reads, 

lKew[q) + 

On the other hand, W{q) is an invariant function: 

Kewiq) h'{q) 
Therefore, the A/'-fold supercharge of the form ()3.1|) is expressed as. 



pm = i-^f n ( ^ + ^(^) - ^^(^) 

(AA-l)/2 9 r ^ \ 



fc=-(Ar-i)/2 
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D. An example 

As an example, we will demonstrate the equivalence between the case P{h) = 2h and 
P{h) = 2h?. In the previous paper js^, we classified them as different cases, namely, case 
(1) for the former and case (3a) for the latter. However, it is easy to see that the latter is 
obtained from the former by the following GL{2, K) transformation: 



h = - . (4.21) 



In this case, the transformation of P and Q defined by Eqs. ()4.3p reads. 



P{hoid) = 2hoid, QiKid) = hhoid + biKid + k 







^ P{hnew) = 2hl^^, Q{Kew) = -^0^^^ - WKew " &2- (4.22) 

From Eq. ()4.4|) . the functions hoid{<i) and hnewiq) are calculated as. 



1 



hoidiq) = q , hnew{q) = — , (4.23) 



q^ 

which is consistent with Eq. (14.211) . The transvectants appeared in Eq. ()4.6|) become, 

H{P) = -3, (P, g)(^) = -4 {SbXid + biKid - bo) 

^ H{P) = -3hl^, (P, Q)W = 4 {bohl^ - hhl^ - 3hhU . (4.24) 

Substituting Eqs. and (jOH) for Eq. (g^l), we have the potentials and as 

functions of h: 



1 , , , .o &2 

1 



V^iKid) = -Kidih^Kid + hiY + ^(6o T ^M)Kid 



8Kid 4 



+ ^{boT^)-R, (4.25a) 
62(60^ SAT) 



+ 4(^0 ^ AT) - P. (4.25b) 

o 4 

Finally, we confirm the invariance of the potentials as functions of q by substituting Eq. ()4.23|1 
for the above 

v^{KM) = v^{K.eM) = Iq' {hq' + 61)' + |(6o t sAf)q' 

bl ± 2A fbo +Af^-1 h 
The functions P(g) and W{q) are calculated as 



+ -^^^ + -iiboT^f) - R. (4.26) 



E{q) = ^ = m = §^ = --, (4.27) 

KM q Kewiq) q 

W{q) = W{q) = -l(b2q^ + biq+-]. (4.28) 
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Thus, the type A TV-fold supercharge of the form ()3.ip reads, 

k=-{M~l)/2 ^ ^ 'i ■i/ 

E. Classification of Models 

For a given P{h), the function h{q) is determined by Eq. (j4.4p and a particular type A 
A^-fold supersymmetric model is obtained by substituting this h{q) for Eq. ()4.fjj) . Since the 
potentials V^{q) and the A/'-fold supercharge are invariant under the GL(2, K) transforma- 
tion, the type A models can be classified according to the inequivalent elliptic integral ()4.4p 
under the transformation. The elliptic integral ()4.4j) can be classified according to the dis- 
tribution of the zeros of P{h), e.g., multiplicity of the zeros. This idea was first introduced 
in Ref. |4^ to classify the sl(2) quasi- solvable models. Under the transformation ()4.3a|) of 
GL{2,M.) or GL(2,C), every quartic polynomial P{h) with real or complex coefficients is 
equivalent to one of the eight or five forms, respectively, as shown in Table Q] In Table IH 



Case 


GL(2,M) 


GL(2,C) 


I 


1/2 


1/2 


II 


2h 


2h 


III 


2i//i2 


2i//i2 


III' 


uih"^ + 1)2/2 




IV 


2ij{h^ - 1) 


2i/(/i2 - 1) 


IV' 


2iy{h^ + 1) 




V 


2/i3 - g^h/2 - 53/2 


2h^ - g2h/2 - 53/2 


v 


u{h^ + l)[{l-k^)h^ + l]/2 





TABLE I: The representatives of P{h) under the GL(2,R) and GL(2,C) transformations. 

^^92,93 £ K according to the transformation group GL{2,K), and u 0, < k < 1, 
g2 — 27 g"^ 7^ 0. In Ref. jj^li^, more general quasi-solvable M-body systems constructed by 
5[(M + 1) generators are classified according to the above scheme and the explicit form of 
the potential for each the cases is shown. The type A models in this article correspond to 
the models for M = 1 (with 6j ±6j for Hjf) in Ref. (i?! . l48l |. So, we do not repeat the 
exhibition of the potentials in this article. 

V. 2-FOLD AND 2ND-DERIVATIVE POLYNOMIAL SUPERSYMMETRY 

All that we have not investigated yet on type A A/'-fold supersymmetry is the anti- 
commutator of Q]^j- and Qj^f, namely, the mother Hamiltonian. For this purpose, it is quite 
instructive to analyze 2-fold supersymmetry. Let us first consider general 2-fold supersym- 
metry, in which 2-fold supercharge is given by, 

P2= p'^ + Wi{q)p + Wo{q). (5.1) 
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In this case, a pair of Hamiltonians = + V^{q) satisfies the intertwining relation 
P2H2 - H+P2 = if and only if, 



w'liq) w[{qf ^ 2Ci 



Wi{q) 2wi(g)2 Wi{qY 
w'liq) w[{qf , 2Ci 



- ^w[{q). 



(5.2a) 
(5.2b) 



_Wi{q) 2wi{qy Wi{qy 

where Ci are arbitrary constants. From the form of V.^{q)^ it is indicated that wl{q) = 
—Wi{q), that is, Wi{q) is pure imaginary in order that V^{q) be real. The above result 
was first reported in Refs. 0] and later reexamined in Ref. |23|. In Refs. |^, the 
anti-commutator of the 2-fold supercharges was given by the following form: 



2H2 = {Ql Q2} = 4 (H2 + C2)' + Ci, 



(5.3) 



that is, it is a polynomial of degree 2 in H2. Later, it was proved in Ref. [20[ that, for all 
J\f, the anti-commutator of Q]^ and Q^v" becomes a polynomial of degree J\f in the original 
Hamiltonian H^v" if Q\f and Qj\f satisfying E g. (12.71) are uniquely determined for the given 
H^v"- Furthermore, it was also proved in Ref. j2y] that, for the above 2- fold supersymmetric 
systems, 2-fold supercharges are unique unless there is a constant C3 which satisfies. 



or equivalently. 



w'liq) - 2tw,iq)w[iq) - 2tV,-'iq) = 2Csw[iq), 



w[{q)-z{w,{q)Y - 2iV,'{q) = 2CsW,{q), 



(5.4) 



(5.5) 



where the integral constant is omitted since it can be absorbed in V2~{q). If Eq. ()5.5|) is 
the case, the 2-fold supersymmetric Hamiltonians satisfy another intertwining relation 
AP2H2 — H2 AP2 = and its conjugation with the following 1-fold supercharges: 



AP2=p + wiiq) - iC3, AP^ =p-wi{q)+ iCs. 



(5.6) 



This result is almost evident from the fact that the condition fj5.5|) together with Eq. ()5.2a|) 
imphes. 



(5.7) 



that is, the 2-fold supersymmetric Hamiltonian H2 in this case is simultaneously ordinary 
supersymmetric (except for the irrelevant constant term). 



[AQ2, H2] = [AQl H2] = 0, {AQl AQ2} = 2H2 + a 
with respect to the supercharges defined by, 

AQ2 = AP^^, AQl = AP2^l 



2 

3' 



(5.8) 



(5.9) 
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Therefore, if we define new 2-fold supercharges by, 

Q2(A) = Q2 + A*Z\g2 = (4 + AMPj)V^, (5.10a) 
QlW =Ql + \AQl= (P2 + A AP2)tpK (5.10b) 

the 2-fold supersymmetric Hamiltonian H2 with ^1(13') satisfying Eq. ()5.5p commutes with 
Q2(A) and QI{X) for arbitrary A: 

[Q2(A),H2] = [g^(A),H2] =0. (5.11) 

On the other hand, the anti-commutator of these new 2-fold supercharges becomes, 

{QlW, Q2{X)} = 4 (H2 + C2)' + |Ap (2H2 + Cl) + Ci 

+ \{Q2,AQI} + \*{qIaQ2}. (5.12) 

It is now evident that the above anti-commutator cannot be, in general, a polynomial in H2 
since it contains the term as (A + X*)p^. Therefore, 2-fold supersymmetry does not always 
correspond to 2nd-derivative polynomial supersymmetry. 

Next, let us return to the case of type A 2-fold supersymmetry. From the definition of 
the type A M-fo\d supercharge ()3.1|) . type A 2-fold supersymmetry is a special case of 2- fold 
supersymmetry where Wi{q) and Wo{q) in Eq. ()5.1|) are given by, 

w,iq) = ~2^W{q), Woiq) = - (w\q) + Wiqf + ^ - ^) . (5.13) 

Substituting the above for the condition ()5.2|) . we see that Eq. ()5.2a|) is equivalent to 
Eq. (IT^ while Eq. (l5^ to Eq. (l3l2E|) . Furthermore, if we substitute Eq. (IKT^ for 
Eq. ()5.5|) . we find that the uniqueness of the type A 2- fold supercharge is guaranteed unless 
there is a constant C3 which satisfies, 

12W{qy - E{qy + 2E'{q) = -IQC^Wiq). (5.14) 

In terms of P{h) and Q{h), the above condition is rewritten as, 

{3Q{hY + H[P]Y = 32ClP{h)Q{hf. (5.15) 

When P{h) is (at most) a polynomial of fourth degree (remember that P{h) can be an arbi- 
trary function in the 2-fold supersymmetric case), both side of Eq. ()5.15|) are (at most) poly- 
nomials of eighth degree belonging to the multiple representation P8, -4 defined in Eq. ()A1|) . 
Therefore, Eq. ()5.15|) is a polynomial identity and the set of its solution constitutes a hy- 
perplane F in the parameter space M^+^ or C^^^ spanned by {aj,6j}. On this hyperplane 
F, 2- fold supercharges are not determined uniquely and the type A 2-fold supercharges can 
be deformed as Eq. (j5.10p without destroying 2-fold supersymmetry. However, it should be 
noted that the polynomiality of the anti-commutator of the undeformed type A 2-fold super- 
charges is preserved on F. The uniqueness of A/'-fold supercharges is only a sufficient but not 
a necessary condition for the polynomiality. Conversely, we can choose A/'-fold supercharges 
such that the anti-commutator of them becomes a polynomial in the Hamiltonian. 

We note that the general 2-fold supersymmetry is weakly quasi-solvable since we cannot 
generally obtain two independent analytic elements of ker P2 where P2 is given by Eq. ()5.1|) 
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with ()5.2b|) . Nevertheless, we can know the two spectra in the solvable sector. The mother 
Hamiltonian ()5.3p is always a polynomial in the original Hamiltonian. Since it corresponds 
to the characteristic polynomial which determines the spectra in the solvable sector, they are 
given by the solutions of 4:{E + C2Y + Ci = 0. This example shows a novel feature of weak 
quasi-solvability. Even if, for a given operator P, there is no analytic element of Vj^f defined 
by Eq. (|2.1ip . we can know the spectra in the solvable sector if the anti-commutator of M- 
fold supercharges constructed from P can be arranged as a polynomial in the Hamiltonian 
H satisfying the weak quasi-solvability condition ()2.12|) . 



VI. MOTHER HAMILTONIANS AND GENERALIZED BENDER-DUNNE 
POLYNOMIALS 

A. Polynomiality of Type A Mother Hamiltonians 

Let us come back to further investigation into type A A^-fold supersymmetry. Suppose the 
pair of type A A/"- fold supersymmetric Hamiltonians = /2 + V^{q) with the potentials 
fl3.2a|) also satisfies the intertwining relation with respect to an A^-fold supercharge (A/" > 
M.) given by. 



M-2 



Pm=P^- tAfW{q)p^'^ + ^ 

n=0 

From a direct calculation, we have, 

[PmHj^ - H^Pm) = (a^W^" - 2^^WW' + 2MVj^' - 2w'j^_^ d- 



(6.1) 



M~2 



2{-t)'-^XW'wk + 2i-^]Vj^^^-''^ + 2Af{' ' \ WV, 



fe=0 

+2 (-ir"[%.y, 

n=k+l ^ ^ 



M-l 
k 



-{M-l-k) 



-{n-k) 



2(-i 



,k-l-M„,J 



gk 



Therefore, Wk{q) {k = 0, . . . , Ai — 2) must satisfy, 
AT. 



w 



M-2 



-W" - W^WW + MVZf 



w 



k-l 



—I 



M-2 

n=fc+l 
■k-\-M 



in+1— fc 



n 



N ' 



-(n-k) 



(6.2) 



(6.3a) 



M 



V, 



-{M-k) 



M - 1 

k 



(A; = l,...,A4-2), 
(6.3b) 



in addition to, 

M-2 



n=l 



+ 2(-i 



^M 



MVj^^^^ + {M- l)^^wvJ 



N 



-{M-l) 



(6.4) 
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The set of the first-order differential equations ()6.3|) can be integrated out in the order from 
Wj\f-2{q) to Wo(g) since all the terms appeared in the r.h.s. of Eqs. ()6.3|) becomes known 
functions depending only on E{q) and W{q) in the order. Thus, we obtain, 

Wk{q) = fk{E{q), W{q)) + Ck {k = 0, . . . , M - 2), (6.5) 

where Ck are integral constants. We can put = without any loss of generality; if 
Cm' 7^ for an tW, we can always split the 7W-fold supercharge as, 

Pm {Cm') = Pm + Cm'Pm', (6.6) 

which means the system is also A^'-fold supersymmetric. Therefore, Wk{q) and hence Pm are 
uniquely determined. On the other hand, if the type A A/'-fold supersymmetric potentials 
fl3.2a|) can be rewritten as the type A A^-fold supersymmetric ones with W{q) replaced by 
J\fW{q)/M., that is, if the following relations, 

V^[Eiq),Wiq)] = V^[Eiq),AfWiq)/Ml (6.7) 

are satisfied, it is evident that the type A A/'-fold supersymmetric system is also type A 
A^-fold supersymmetric with respect to the following type A A^-fold supercharge: 

{M-l)/2 . X 

Pm= n [p-tj^W{q)+tkE{q)]. (6.8) 

k=-{M-l)/2 ^ ^ 

As we have just shown, however, Pj^ is unique and thus the A^-fold supercharge determined 
by the solutions of Eqs. (|6.3p must be the type A A^-fold supercharge Eq. (j6.8p . Furthermore, 
the constraint ()6.4|1 must be equivalent to the relation ()6.7|1 . Summarizing the result, we 
have shown that type A A/'-fold supercharge is unique unless the system satisfies Eq. ()6.7p 
for a.n M. (0 < A4 < A/"). The relation ()6.7|) is equivalent to, 

12W{qf - {E{qf - 2E'{q)) = {2Q{hf + M'H[P]) = 0, (6.9) 

and its solutions again constitute a hyperplane T in the parameter space spanned by {oj, bi}. 
Outside the hyperplane, the A/'-fold supercharge is unique and the mother Hamiltonian is 
expressed as a polynomial V_\f in the original Hamiltonian H^^^: 

'HM=l{Ql,QM}='PM{'iiM) for {a„6,}eR5+3\rorC^+='\r. (6.10) 

On the other hand, as has been shown in Eq. ()2.19j) . this polynomial Vj\f corresponds to 
the characteristic polynomial which determines the spectra in the solvable sector |22|. Since 
the characteristic polynomial itself is constructed from a finite algebraic operation, it has 
no discontinuity in the parameter space. Therefore, Eq. ()6.1()|1 must be held on F and thus 
the type A mother Hamiltonian must be a polynomial in the original Hamiltonian H^v" in 
the whole parameter space. 
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B. Generalized Bender Dunne polynomials 



In 1996, Bender and Dunne introduced a set of polynomials which determines the spectra 
in the solvable sector of the quasi-exactly solvable model of case II, namely, the sextic 
anharmonic oscillator js^. Soon after, Finkel et al. generalized the idea to all the sl{2) 
quasi-exactly solvable models [s^. In the following, we will make a generalization without 
imposing the normalizability of the solvable sector. Let the bases {0^} of the solvable 
subspaces be. 



h"^-' {G^U^y' (n=l,...,Ar), 



and the gauged Hamiltonians be, 



dh 



Then we have. 



-{n - l){n - 2)P{h)(p. 



± 

n-2 



From Eqs. ()3.41|) . and ()3.50b|) with ()3.52p and ()3.55p . the P^{h) are expressed as, 

J\f-2 



pm 



-P'{h)TQ{h), 



(6.11) 

(6.12) 
(6.13) 

(6.14a) 
(6.14b) 



Substituting Eqs. ()3.42|) . ()3.43|) and ()6.14|) for Eq. ()6.13|) . we obtain the matrix elements 
defined in Eq. dTTl) : 



n,m 



-{n -U){n-U+ l)a4 0^+2 - {n - U) 



{n-l){n-U) + UM-l){N-2) 
6 



n - y ) ag =F 62 



«2 =F I — \ bi + R 



{n-l) 



M + 2 



n 



ai =F 60 



_i - (n - l)(n - 2)ao 



^2- 



(6.15) 



As we have discussed in SectionHV| there are five independent type A A/'-fold supersymmetric 
models under the GL{2, C) transformation. By a suitable GL{2, C) transformation, we can 
always transform P{h) so that = in all the five cases, see Table ITTl In Table UTl Cj {i = 
1, 2, 3) denote the three different single roots of the algebraic equation 2x^ —g2x/2 — g3/2 = 
and hi {i = 2, 3) are given by hi = Cj — ei. So, we set ao = hereafter. We note that case 
I corresponds to the case where P{h) has a quadruple root and thus P{h) = constant 
(quadruple root at infinity) or P{h) oc {h — h^Y (quadruple root at finite ho). Therefore, 
in contrast to in all the other cases, case II-V, we cannot simultaneously set 04 = and 
ao = by any GL{2, C) transformation in case I. We also drop another irrelevant constant 
by putting R = 0. 
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Case 


h 


P{h) 


I 


l/h 




II 


h 


2h 


III 


h 




IV 


{h + 2)/h 


2uh'^{h + l) 


V 


{eih + h2h^)/h 


2h(h - h2){h - hs) 



TABLE II: The GL(2, C) transformations h which convert the standard forms of P{h) in Tabled 
into other forms P{h) satisfying oq = 0. 



In the next, we introduce a set of functions Pii^{E) by putting a solution of the 
Schrodinger equation as follows: 



^(±6o)"n! 



n=0 



0n+l 



-ai)" n! r I n 



AA-2 bo 
2 "T ai 



(6.16) 



From the requirement that the above ip"^ satisfies H^ij)^ 



if ai = 0, 
if ai ^ 0. 

Eip^, we obtain a four-term 



recursion relation for P]t^{E): 

Pi:^AE) = {E + A\^)Pi^\E) - n{n - Ar)5[^Pi^,(E) 



where A\^, and Cn^^ are given by. 



n(n-X+l) + l{Af-l){^- 2) 
6 

n - — ] aiTbo 



0-2 =F - 



AT 



- TT Ms =F 02 



a4 



n - — ] aiTbo 



n 



ai =F &o 



(6.17) 

(6.18a) 
(6.18b) 
(6.18c) 



We can set P^\e) = 1 without any loss of generality. Then, each pIt^\e) generated by 
Eq. ()6.17|) becomes a monic polynomial of nth degree. We call it a generalized Bender- 
Dunne polynomial (GBDP). In the case of = 0, Eq. ()6.17|) reduces to a three-term 
recursion relation with (in general) n{n — Af)Bl{^ ^ for n ^ Af and thus the set of 
p}i^{E) forms a weakly orthogonal polynomial system ji^. Therefore, case I is special in 
the sense that the set of polynomials Pjt^ does not form a weakly orthogonal polynomial 
system. One of the most notable properties of the GBDP is the factorization property; due 
to the structure of Eq. ()6.17|) . all the polynomials P]t^{E) for n > Af are factorized as, 

Pi^liE) = Q^J^Ke)p!^'\e) (n>0), (6.19) 
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where Q\^{E) is a polynomial of degree n satisfying another four-term recursion relation: 

+ n{n -l){n + Af){n + Af- l)C^„gL^U^)- (6-20) 

In the case of 04 = 0, Eq. ()(i.2()j) reduces to a three-term recursion relation with (in general) 
n{n+Af)B^J^^ 7^ for all n > and thus the set of Q\^{E) forms an orthogonal polynomial 

system j4^. Therefore, case I is again special in the sense that the set of polynomials 
does not form an orthogonal polynomial system. Due to the factorization property ()6.19p . 
p)^\e) reserves special status among the GBDPs. We call P^{E) Af-th critical GBDP 
after the terminology in Ref. [s^ . 

When E takes one of the spectral values En {n = 1, . . . ,J\f) in the solvable sector, 
must be an element of V^. From the factorization property ()(119j) . the condition ijj^ G 
is fulfilled if and only if, 

pjf{En) = {n = l,...,Af). (6.21) 

This means that all the zeros of the A/'-th critical GBDP must correspond to the spectra in 
the solvable sector of the quasi-solvable Hamiltonians H^. On the other hand, they are also 
given by solutions of the characteristic equation ()2.15|) . Therefore, each of the critical GBDP 
is proportional to the corresponding characteristic polynomial. Comparing the coefficients 
of the highest-degree term with each other, we have, 

det M%{E) = 2^Pi^(E). (6.22) 

As we have proved before, the type A mother Hamiltonian is expressed solely by the char- 
acteristic polynomial of the original Hamiltonian as Eq. ()2.19p . Combining Eq. ()2.19|) with 
()6.22|) . we finally obtain an intriguing relation: 

nM=^ det M± (Hat) = 2^-ipi?^(H^). (6.23) 

Furthermore, remembering that the mother Hamiltonian is defined by the anti-commutator 
of the A/'-fold supercharges (j2.16|) . we get the complete type A A/'-fold superalgebra: 

[Qm,^m] = [Q]^,H^] =0, (6.24a) 
{Qn, Q^} = {Q]v^, QJv-} = 0, (6.24b) 
{Q]^,QAr}=2-^P^](H^). (6.24c) 

C. Examples 

In order to confirm the previous argument, we will show the explicit results for 
Af = 1, 2, 3. By solving the recursion relation ()6.17|) . we obtain pIt^\e) {n < Af) as follows: 

1) Ar= 1: 

pI^\e) = E, (6.25) 
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2) AT = 2: 



2 ' 



ft2 



3) AT = 3: 
pf\E) 
pf\E) 

pf\E) 



E+^^±h, 
E^-^{a2TSh)E 



(4a2 — 9aia3 =]= 180162 ± 1202^1 =F ISfls^o ^ 366062) , 

18 

+ - (3aia3 -al + I26062 - 36i) E - — [2al - gaiasag 
3 27 



(6.26a) 
(6.26b) 

(6.27a) 

(6.27b) 



+270^04 - 1080460 + 54036061 - 180261 - 36026062 + 54016162) . (6.27c) 
On the other hand, the direct calculation of the mother Hamiltonians reads as follows: 
1) Ar= 1: 

27^1 = 2Hi, 



2) A/" = 2: 



3) A/" = 3: 



27^2=4(H2)' + /^2[Q], 



16 



27^3 = 8 (H3)' - -{t2[P] - 3A[g])H3 + -{UP] + 9/i,2[P, Q]). 



(6.28) 



(6.29) 



(6.30) 



In the above, D2, ^2, Js and /i^2 are the absolute invariants expressed in terms of E{q) 
and W{q), see Eqs. flMil) - (p^ . Comparing the critical GBDPs, Eqs. (jFT^ . (lOHE]) . 
and (j6.27c|) . obtained by solving the recursion relation (|6.17p with the mother Hamiltonians 
fl6.28|) - (j6.30|) . and noting that we have put oq = 0, we confirm the relation ()6.23|) for Af = 
1,2,3. We also note that since the mother Hamiltonians are invariant under the GL{2,K) 
transformation, all the coefficients of the critical GBDPs should be expressed solely in terms 
of the absolute invariants listed in Eq. ()A6|1 . as the above examples indicate. From the 
facts that all the GBDPs are homogeneous polynomials in E, ai and 6j due to the structure 
of the recursion relation (j6.17|) . and that all the critical GBDPs are symmetric under the 
transformation hi —hi, the general form of the critical GBDPs should be. 



P^^{E) = E'^ 



k=0 



^2ni+2n2+3n3+3n4+4n5+12n6,Ar-fcX 



.ni,...,n6 



X C, 



fc;ni,...,n6 K^'i) 



(J3)"^ (/l,2)"^ (/2,2)"-^' (/3,3)'"^ 



(6.31) 
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where C^^^ „g are constants. For example, the first five are, 

p'\^\e) = E, (6.32a) 
pf\E) = E^ + \D,, (6.32b) 

Pf^ (E) = E'-^it2-3D2)E + ^ (j3 + 9/1,2) , (6.32c) 

Pj'l (E) = E^-l- (422 - 5D2) E^ + 4/1,2 E+{t2f + 2/2,2 + t2D2 + ^{D2)\ (6.32d) 
/ ib 

Pf 1(E) = E'- (7z2 - 5/^2) - 2 (j3 - 7/1,2) E^ 

+ 4 (3 (22)' + 4/2,2) E + 8it2- D2) (j3 + h,2) ■ (6.32e) 

Since the zeros of the critical GBDPs correspond to the spectra of in the solvable sector, 
these spectra can be regarded as functions of the six absolute invariants. Furthermore, as 
an interpretation of Eq. ()6.31|) . the critical GBDPs, and equivalently the type A mother 
Hamiltonians, can be regarded as generating functions of the absolute invariants composed 
of P{h) and Q{h). 



VII. CONCLUDING REMARKS 

In this article, we have fully investigated general aspects of type A A/'-fold supersym- 
metry. The two different approaches in Section IIIII reveal both the analytic and algebraic 
structures of the systems. The intimate relation between the algebraic forms of the pair 
Hamiltonians, Eqs. (j3.38p and (j3.57|) . provides an interesting problem. Suppose we have a 
quasi-solvable system H which is represented by a quadratic polynomial of a set of first-order 
differential operators constituting a finite dimensional representation of a Lie algebra. If we 
construct another system H from H by interchanging the order of the quadratic terms and 
interchanging the sign of the coefficients of the linear terms, do H and H (with suitable 
gauge transformations) always form an A/'-fold supersymmetric pair? This problem may be 
extended to quasi- solvable many-body systems j47|, liSf • 

The invariance of the Hamiltonians as well as the type A jV-fold supercharge under 
the GL{2,K) transformation enable us to obtain the complete classification of the type 
A models. This invariance also plays an essential role in generalizing the Bender-Dunne 
polynomials to for all the type A models. With regard to the issue, we should stress that, up 
to now, we have found no reason that a set of the polynomials should be weakly orthogonal, 
although all the existing papers concerning about the issue HHHHHlHlIll, as far 



as we know, have required weak orthogonality by imposing, for instance, the normalizabiliy, 
or restricted the considerations in which weak orthogonality is fulfilled. One of the most 
essential properties that sets of polynomials associated with quasi-solvable systems should 
share is the factorization property ()6.19|) . The most general form of a recursion relation for 

a set of polynomials Pjf^ which guarantee the factorization property may be the following: 



K 



P'fAiE) = {E + A^^^)Pi^\E) + Y: 



k=l 



/Yiin-l)in-Af-l)A^^jpniE) 
1=0 



(7.1) 



Then, an interesting question is, what conditions a set of the coefficients {^4^} should satisfy 
for the existence of a quasi-solvable system whose spectra in its solvable sector are given by 
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the zeros of the critical polynomial (E) obtained by the above recursion relation. This 
question may provide an alternative way to find out a new quasi-solvable and also a new 
A/'-fold supersymmetric model. 
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APPENDIX A: INVARIANTS OF A SYSTEM OF POLYNOMIALS 

Since the type A TV-fold supersymmetric models have the underlying GL{2, K) invariance 
discussed in Section IIVI all the relevant quantities should be expressed in terms of the 
covariants and invariants under the transformation. In this appendix, we summarize the 
covariants and invariants of a system of polynomials needed in this article. For more details, 
see Ref. ji^ and references cited therein. 

The multiplier representation p^.i of GL{2,K) on the space of polynomials of degree at 
most m is defined by, 

F{h) ^ F{h) = A\^h + 6rF{h), (Al) 

where h and A are defined by Eq. ()4.H) . For example, we see from Eqs. ()4.3p that P{h) 
belongs to the representation P4, -2 while Q{h) belongs to the representation p2,-i- 

In order to construct a complete system of covariants and invariants, the process of 
transvection is essential. Let F{h) be a polynomial of degree at most m belonging to 
the representation j and G{h) be a polynomial of degree at most n belonging to the 
representation Pn,j- The rth transvectant of F and G is defined by. 



(A2) 



Then, (P, G)^^^ is a polynomial of degree at most m + n — 2r belonging to the representation 

Pm+n—2r, i+j+r- 

For a single quadratic polynomial Q, there are two independent covariants, namely, Q 
itself and its discriminant D2 given by, 

D2[Q] = liQ,Qy'\ (A3) 

For a single quartic polynomial P, there are five independent covariants, namely, P itself, 
the Hessian H and the Jacobian J given by, 

H[P] = ^(P, Pf^ = PP" - I {P'f , (A4a) 

J[P] = ^(P, P-)(1) = J4p2p(3) _ Qpp'p" + 3 (p')3j ^ (^4^^ 
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and the two invariants given by, 

^2[P] = MP] = ^(p^hY^^- (A5) 

The discriminant of P is expressed in terms of the above invariants as -Del-P] = '^2[-P]^~j3[-P]^- 
According to the invariant theory of polynomials, the complete list of independent 
covariants for the system consisting of a quartic and a quadratic polynomial is the following: 

1) Absolute invariants po.o^ 

D2[Q], ^2[P], J3[P], /l,2[P,Q], /2,2[P,Q], /3,3[P,Q], (A6) 

2) Quadratic covariants P2,-i- 

Q, {P,Qf\ {H,QY'\ {P,QT\ {H,QT\ iJ,QT\ (A7) 

3) Quartic covariants P4, _2: 

p, H[P], (p,g)«, {h,qY'\ {j,qY'\ (A8) 

4) Sextic covariant pe, -3: 

J[P]. (A9) 
In the above, the absolute invariants /i,2, -^2,2 and /s^s are given by, 

h,2[P,Q] = ^{P,QT\ h2[P.Q] = ^iH,QT\ ^3,3[P,Q] = ^(J,g')('^ (AlO) 

We can express the above quantities in terms of h{q), E{q) and W{q) with the aid of the 
following relations derived from Eqs. ()3.3U|) and ()3.44|) : 

P = ^{hf, P' = h" = Eh', P" = E' + E'', (Alia) 
p(3) = 1 (E" + 2EE') , P(^) = [^(3) + EE" + 2 (E'? - 2E^E'] , (Allb) 



Q = -Wh\ Q' = - {W + EW) , (A12a) 
Q" = -1 (W" + EW + E'W) . (Al2b) 

For example, the Hessian H and the Jacobian J are expressed as, 

H[P] = ^ {2E' - E') , J[P] = {E" - EE') . (A13) 

Among the covariants ()A6|) - ()A9|) . the absolute invariants play an essential role in Section IVll 
In the following, we show the explicit forms of them in terms of the coefficients of the 
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polynomials P{h) and Q{h), as well as in terms of E{q) and W{q) (for the first four invariants 
because the expressions for the last two are lengthy): 



D2[Q] = 2QQ" - {Q'f 

= 2WW" - {W'f + {2E' - E^) 
= 46062 - bl 



(A14) 



2pp(4) _ 2P'p(3) + (^p"y 
^(3) _ EE" + 3 {E'f - AE'^E' + E' 



12aoa4 — 3aia3 + ag, 



(A15) 



2j3[P] = ^ 12PP"P(^) - 6P (P(3))' - 9 (P')' P(^) + 6P'P"P(3) - 2 (P")' 
8 

1 



8 



6E'£;(3) _ 3 (^E" f - SE^E^^^ 

+10 + 3E^E" - 2AE^ {E'f + 12£;^£;' - 2E^ 
720002^4 — 21aQa\ — 270^04 + 9010203 — 2o2, 



(A16) 



47i,2[P, Q\ = P^^^Q' - 2P(3)QQ' + 2P" [gg" + (Q')'] - 6P'Q'Q" + 6P (Q")" 
^(3) _ EE" + 7 {E'f - SE'^E' + 2E^] W"^ 

- 2 {E" - EE') WW' + {2E' - E'^) UwW" + {W'f \ + 3 {W"f 
— 4 (60463 — 3036061 + 20260^2 + Q2&1 — 3016162 + 60062) ) (Al''') 

2/2,2[p, Q] = \ {i/(^)g' - 2i/(3)gg' + 2h" \qq" + (g')'] - ei/'g'g" + &h (g")'} 

= 3 (802O4 — 303) 60 ~ 6 (601O4 — 0203) hobi 

+ (2400O4 + 30103 — 202) (26062 + ^1) ~ 6 (600O3 — O1O2) 6162 

+ 3 (800O2 - 3oi) 62, (A18) 

/3,3[P, Q] = l^ {^^'^Q' - 3J(')g'g' + 3j(^)g [gg" + 2 (g')'] - 6j(=^)g' [3gg" + {Q'f 
+isj"Q" [gg" + 2 (g')'] - Qo/g' (g")' + 90J (g")'} 
= (80104 — 4020304 + 03) 60 — (160004 + 2010304 — 40204 + 0203) 6061 

+ (800O3O4 — 401O2O4 + 01O3) (6062 + 606^) — 6 (00O3 — 0^04) 606162 

- (00O3 - 0^04) b\ — (800O1O4 - 400O2O3 + 0^03) (6062 + 6^62) 

+ (I600O4 + 200O1O3 — 400O2 + 0^02) 6162 — (800O3 — 4ooOi02 + o^) 62- (A19) 
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